Casimir-Polder forces on moving atoms 
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Polarisable atoms and molecules experience the Casimir-Polder force near magnetoelectric bod- 
ies, a force that is induced by quantum fluctuations of the electromagnetic field and the matter. 
Atoms and molecules in relative motion to a magnetoelectric surface experience an additional, 
velocity-dependent force. We present a full quantum-mechanical treatment of this force and iden- 
tify a generalised Doppler effect, the time delay between photon emission and reabsorption, and the 
Rontgen interaction as its three sources. For ground-state atoms, the force is very small and always 
decelerating, hence commonly known as quantum friction. For atom and molecules in electronically 
excited states, on the contrary, both decelerating and accelerating forces can occur depending on 
the magnitude of the atomic transition frequency relative to the surface plasmon frequency. 

PACS numbers: 12.20.-m, 34.35. +a, 42.50.Nn, 42.50.Wk 
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I. INTRODUCTION 

The ground-state fluctuations of the electromagnetic 
field lead to several inherently quantum effects such as 
the spontaneous decay of excited atoms and molecules 
as well as dispersion forces Forces between isolated 
atoms that are mediated by the quantum vacuum are 
known as van der Waals forces while forces between 
macroscopic bodies are referred to as Casimir forces ||- 
The third type of dispersion (in a sense an interpolation 
between these two extreme cases) is the Casimir-Polder 
(CP) force exerted on single atoms near macroscopic bod- 
ies [|. 

For a two-level atom with transition frequency uj a and 
electric dipole moment d located at a distance z A away 
from a perfectly conducting plate, the short-distance 
(non-retarded) Q and long-distance (retarded) || lim- 
its of the CP potential take the well-known forms 



U mc t(z A ) 
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These potentials, acting on atoms at rest, lead to con- 
servative forces perpendicular to the plate's surface. CP 
forces (as well as all other dispersion forces) play im- 
portant roles as limiting factors in efforts to miniaturise 
atom-optical devices [6| , and have been measured at dis- 
tances as small as 6/im Q- 

Casimir-Polder forces are well understood far beyond 
the aforementioned special case of a perfectly conducting 
date, with magnetoelectric materials of arbitrary shape 
8| and finite temperature being investigated theoretically 
9f.[i"c| as well as experimentally [ll|]. While most theoret- 
ical investigations are based upon Lifshitz' macroscopic 
treatment [ill or a linear-response description [3], full 
quantum theories based upon electromagnetic-field quan- 



tisation in magnetoelectrics have also been developed [8|] . 
In the latter approach, the operator-valued Lorcntz force 



d 3 r 



P A (t)&(t)+3 A (t) x B(r) 



(2) 



on the atomic charge and current densities due to the 
body-assisted electromagnetic fields is computed. In 
the long-wavelength approximation, it leads to the well- 
known expression 



d-E(r A ) 



d 
dt 



d x B(r A ) 



(3) 
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with d denoting the atomic electric dipole moment oper- 
ator. 

Intuitively, dispersion forces can be understood as 
dipole-dipole forces generated by spontaneous polarisa- 
tion due to the electric-field fluctuations. Therefore, CP 
forces on atoms at rest act either towards (attractively) 
or away (repulsively) from the macroscopic body. For 
atoms in motion, retardation delays the dipole response, 
and a force component emerges along the direction of 
motion. In most cases, this force component acts against 
the motion and is thus the origin of quantum friction. 

Quantum friction forces have traditionally been stud- 
ied within a linear-response formalism [3, [l5|, [l6|, [l7l |. 
Evaluating the correlated quantum fluctuations of mov- 
ing atom and dielectrics, the friction force on ground- 
state atoms can be obtained in this way. However, the 
predicted forces are typically very small. For the station- 
ary case it is well known that CP forces can be resonantly 
enhanced for excited atoms [1, [HI . For such nonequilib- 
rium situations, linear-response methods cannot be ap- 
plied and a more detailed investigation of the atom-field 
dynamics becomes necessary. 

In this article, we develop a full quantum theory of 
the velocity-dependent CP force. In particular, we will 
show that for atoms and molecules in electronically ex- 
cited states, both decelerating and accelerating forces can 
occur depending on relative magnitude of the frequencies 
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of the atomic transition and the surface plasmon. The 
article is organised as follows: After briefly reviewing the 
formalism of macroscopic quantum electrodynamics in 
Sec. UH we study the atom-field dynamics in Sec. IIII Al 
before investigating the resulting force in Sec. IIII Bl and 
applying our results to the quantum friction scenario in 
Sec. IIII CI We illustrate the theory with representative 
examples in Sec. IIVI followed by a summary in Sec. [V] 



II. BASIC FORMULAS 

Let us assume an arbitrary arrangement of dispers- 
ing and absorbing magnetoelectric bodies, characterised 
by their complex- valued, Kramers-Kronig consistent per- 
mittivity s(r, u>) and permeability /i(r, u>). The Hamilto- 
nian of the quantum electromagnetic field and the bodies 
can be given as (for a recent review, see Ref. [l9| ) 



H F = 



E 

A— e.m 



d 3 7 



duhuji{{r,u}) ■ fx(r,cj) (4) 



where the fundamental variables f|(r, and f A (r, u>) are 
creation and annihilation operators for the elementary 
electric (A = e) and magnetic (A = m) excitations of the 
system; they obey the bosonic commutation relations 

f A (r,^),fj,(r>')] =S X yd(r-r')S(cj-u J '). (5) 

The electric and magnetic fields can be expanded in terms 
of the fundamental variables according to 



E(r) 



dwE(r,w) + H.c 

o 

E 

A=e,m 

+ H.c 



d 3 -' 



;- / duG A (r,r» • f A (r',u;) 
o 



(0) 



B(r) = / dwB(r,w) + H.c. 
Jo 

J dV l°° x G x (r,r',u>)-h(r',iv) 

(7) 



A— e,m 

+ H.c. 



with coefficients 



C e (r,r' )U ) = i^/llm E (r»G(r,r», (8) 



G m (r,r',w) = i-W j-r-, — rry V'x G r', r, w) . 

(9) 

Here, G is the classical Green tensor as uniquely defined 
by the inhomogeneous Hclmholtz equation 



V x — V x - — e(r, w) 



G(r,r',w) = <5(r-r') 
(10) 



together with the boundary condition 

G(r,r\u)->0 for |r - r'| ->• oo. (11) 

The Green tensor is an analytic function in the upper half 
of the complex frequency plane and it has the following 
useful properties: 

G(r,r',-u ) *) = G*(r,r',oj), (12) 
G(r',r,c) = G T (r,r',c), (13) 

Y, J d 3 S G A (r,s,^).G A T (r',s,c) (14) 



A— e.m 



cu ImG(r,r',w) 



(15) 



The Hamiltonian describing the internal dynamics of 
an atom with eigenenergies E n and eigenstates \n) can 
be given as 



Ha — Y^ EkAkk 



(16) 



(Akk = \k)(k\: atomic flip operators). Throughout this 
article, we will assume that the centre-of-mass motion is 
sufficiently slow so that it separates from the internal dy- 
namics in the spirit of a Born-Oppcnhcimcr approxima- 
tion. The interaction of the atom with the body-assisted 
electromagnetic field is then adequately described by the 
atom-field coupling Hamiltonian for given centre-of mass 
position ta and velocity v^, which in the multipolar cou- 
pling scheme and electric-dipole approximation reads 

H AF = - d • E(r A ) - d • w A x B(r A ) 

= - y^dfc; ■ E(r A )A k i - Y d kl ■ va x B(r A )A kl . 



k-i 



(17) 



The first term is the familiar electric-dipole interaction 
while the second term is the Rontgen interaction as- 
sociated with the centre-of-mass motion. Combining 
Eqs. (g]), fig]), and (JT7J), the total Hamiltonian of the 
atom-body-field system reads 



H = H A + H F + H 



AF- 



(18) 



Finally, the total Lorentz force on the atomic charge 
and current distribution can in electric-dipole approxi- 
mation be given as 



F = V A [d • E(r A ) + d • v A x B(r A )] + — [d x B(r A )] 
= V A Y t dfci ' M t a)^u + d M ■ x B(r A )A kl ] 



+ di Et dfei x B(r A )ifei]. 



(19) 
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III. CASIMIR-POLDER FORCE ON A MOVING 
ATOM 

The Casimir-Polder force on an atom is the quantum 
average of the Lorentz force (fH)|) with the body-assisted 
field being in its ground state. To evaluate this expres- 
sion, we first need to solve the coupled atom-field dy- 
namics. 



A. Atom— field dynamics 

Using the Hamiltonian (|18|) . the Heisenberg equations 
of motion of the atomic and field operators are found to 
be 



A mn — iWnmAiJB + y ] (d n fc-^-mfc ~~ ^km^kn) ' ^(f^) 
k 

+ j ^2(d n kA mk - d km A kn ) -v A xB(r A ) (20) 



and 



f A (r,w) = -iwf A (r,w) + ^^G A T (r,rA,w) ■ d kt A k i 

k,l 

+ ^]T[G A T (r,r A , W ) x V'] x v A ■ d kl A kl (21) 

k,l 

(by convention, V and V only act on the first or second 
argument of the Green tensor, respectively). The latter 
equation is formally solved by 

fx(r,w,t) = fA,f(r,w,t)+f Aia (r,w,t) (22) 

where 

& jf (r,w,t) =e- iwt f A (r, W ), 



(23) 



f AjS (r,w,f) 



^ / t dre- i -^G A T [r,r A (r),c] ■ d kl A kl (r) 

k,l 

- / dre^^Cf [r,r A (r),w] x V'} 



x • d fci ^4 fe i(r) 



(24) 



determine the free and source parts of the electromag- 
netic field. 

We assume that the atom moves with uniform, nonrel- 
ativistic speed (v A "C c) and we are seeking a solution to 
the system of Eqs. (|20|) and (f2Tj) within linear order of 
Va- We may hence write 



i-a(t) = r A (t) -(t- t)v a ; 



(25) 



and after substituting Eqs. ([22 ]) -([2"3] ) into Eq. ((6|), using 
the integral relation (|14p . and applying a linear Taylor 



expansion in v_a, the time-dependent frequency compo- 
nents of the electric field are given by 



E(r,w,f)=E f (r,w,t)+E B (r,w,t) 



with 



E f (r,w,t) = e- iwt &(r,co) 



(26) 



(27) 



^ W 2 f dr e"^ 4 -^ V Im G(r, r A , w) ■ d kl A kl (r) 

11 J o tr 



^u, 2 / drft-rje- 1 "^ 
n Jo 



x 



^ImG(r, r A ,a;)-d w (V'-v A )i w (r) 



k.l 



^ w / dr e - iw (*- T) V Im [G(r, r A , w) x V'l 



x v A ■ d k iA k i(r) 



(28) 



[r A = r A (t)]. The magnetic field ([7]) only enters the 
equations of motion in conjunction with a factor v^, so 
we only require its zero-order expansion in the velocity: 

B(r, u, t) = Bf (r, w, i) + B s (r, w, i) (29) 

with 

Bf(r,w,t) =e- iwt B(r,w) (30) 
B s (r,w,i) 

^VxImG(r,r A ,w)-d fc |4(T). (31) 



Mo 

7T 

X 



t 

/ dre^*^ 
o 



fe.i 



We can next substitute our solutions (|27|) - ([3lj) for the 
time-independent electromagnetic fields into the equa- 
tion of motion (|2"U)) for the atomic flip operators. Noting 
that the total field operators E(r, u>, t) and B(r, u, t) com- 
mute with the atomic flip operators at equal times, we 
arrange all products such that creation operators f| (r, oj) 
are always at the left and annihilation operators f A (r, uj) 
are always at the right. Assuming the field to be initially 
prepared in its vacuum state and taking expectation val- 
ues, all contributions from the source fields vanish. For 
weak atom-field coupling, the time integrals can be eval- 
uated with the aid of the Markov approximation, 



dTe-^ t - T \A l0 (t)A kl { T )) 



{A %3 {t)A kl {i)) 
= (A u (t))5 jk 



dTe -i ("~ LD " i ^*~ r ) 
V 

it5(uj - uik) - i — 

LU - UJl k 



(32) 
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(V: principal value); similarly we have 



dr (t-T)e- i -( t -^(4-(t)i fc; (^)> 



(A u (t))S jk 



V 



UJ - LU[ k 



+ mS(ui - Cbik) 



where the shifted atomic transition frequencies 

&mn = Umn + Su! m — SuJ„ 



(33) 



(34) 



have yet to be determined. 

For a nondegenerate atom, the resulting equations of 
motion for the off-diagonal atomic flip operators decou- 
ple from each other as well as from the diagonal ones. 
In addition, we consider an atom whose internal Hamil- 
tonian (|16p is time-reversal invariant, so that we may 
assume real dipole-matrix elements. After a lengthy, but 
straightforward calculation, we finally obtain the follow- 
ing equations of motion for the internal atomic density 
matrix elements a mn — (A nm ): 



Pn 



kPk, 



(35) 



ow = [-i&mn - (r m + T n )/2]a mn for m ^ n, (36) 

where we have introduced the probabilities p n = o nn . 
The equations of motion for a moving atom have exactly 
the same form as for an atom at rest: The population 
of the diagonal density matrix elements is governed by 
spontaneous decay, while the off-diagonal ones undergo 
damped oscillations. However, the respective transition 
rates 



(37) 



(the primes indicate derivatives with respect to w). Here, 
we have decomposed the Green tensor into its bulk (free- 
space) and scattering parts according to 

G{r,r',oj) = G ( -°\r,r',u>)+G^(r,r',u ) ) (45) 

and have discarded the Lamb-shift contribution due to 
G'- -' from the frequency shift Suj k (r A ) (as the free-space 
Lamb shift is assumed to be already included in the bare 
transition frequencies uj mn ). We have further exploited 
the symmetry fp~3|) of the Green tensor which implies 



^Va, showing that the translationally invariant 
bulk Green tensor does not contribute to the velocity- 
dependent shifts and rates. 

It is worth noting that for real dipole matrix elements, 
the contributions from the Rontgen interaction exactly 
cancel. As a result, the velocity dependence of these 
quantities is entirely due to the fact that the moving atom 
emits and receives the electromagnetic field at different 
positions. The velocity-dependent contributions are pro- 
portional to the total derivative of the scattering Green 
tensor along the direction of motion. As a consequence, 
the decay rates and frequency shifts are unaffected by 
uniform motion in a direction along which the environ- 
ment is translationally invariant (e.g. motion parallel to 
a plate or a cylinder). It is worth pointing out that such 
a vanishing of velocity-dependent frequency shifts does 
not necessarily imply that the velocity-dependent part of 
the CP force must also be zero. One should bear in mind 
that all of the above has only been shown within linear 
order in the velocity. 



Casimir— Polder force 



and frequency shifts 



Sui n = ^2 SuJ n 



are affected by the atomic motion: 
Mo 



irn 







U nk 



Mo 



6w*(r A , v A ) = §r e(£ nfc )(vA • V A ) 

x [w 2 d„ fc • lmG(r A ,r A ,uj) ■ dkn] u=S)t 

r* = r*(r A ) + r*(r A ,v A ), 



r£(rA,VA) = -^(vA- V A )V /°°d w 

Wd nk ■ lmG w (r A ,r A ,uj) ■ d kn Y 
x — 

Writ — W 



(38) 



(39) 

Lu 2 d nk ■ ImG (1) (rA,rA, w) ■ d kn 



(40) 



(41) 
(42) 



F^(i-a) = ^ @{ij nk )u>l k d nk ■ ImG(rA,r A ,<I>„ fc ) • d 



fen ) 

(43) 



(44) 



Having solved the coupled atom-field dynamics, we 
can now evaluate the quantum average of the Lorentz 
force (fH?|). We restrict our attention to the pure disper- 
sion force by again assuming the field to be initially pre- 
pared in its ground state. The atom may initially be in 
an arbitrary incoherent superposition of internal energy 
eigenstates. For an atom at rest, it has been found that 
the third term in Eq. (|19p . which involves a total time 
derivative, does not contribute to the force on atoms in 
incoherent internal states. We have explicitly checked 
that the same is true here, so that we only need to con- 
sider the force 

F = V A (d ■ E(fa) + d ■ va x B(r A )> 

du ■ va x (B(r A )A k i)] . 

(46) 



/,■/ 



We begin by substituting the time-dependent electro- 
magnetic fields (|26l) - (f3Tj) where again we retain only 
terms up to linear order in va and we arrange all prod- 
ucts such that the contributions from the free fields van- 
ish. The source fields give rise to intra-atomic correlation 
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functions. By virtue of the quantum regression theorem, and similarly 
Eq. p6p implies that the relevant correlation functions 
are of the form 



(A nk (t)A ln (r)) = 4ie in -(*- T V n „(r) (47) 



with 



(t)S k i 



Qnk = w„fe + i(r„ + r fc )/2. 



(48) 



(lo - Q nk ) 



(50) 



approximation: 



We evaluate time integrals in the spirit of the Markov Again assuming real dipole matrix elements, the resulting 

expression for the CP force can be written in the form 

F(<)=J>„(t)F B ( 51 ) 



dTe -™(t-T)(A nk (t)A ln (T)) 



<J n n(t)hl / dre 



-i(w-n„ k )(t-r) _ a nn(t)iSkl 

ca — f2 n fc 

(*) with 



Vd» fc ■ Im (r A , r^, a;) ■ d fctI 

LO — Q nk 



V(v A ■ VQdnfc ■ lmGW(rA,rA,ij) ■ d fc » 



ift^fj (V - V)d„ fe • v A x [VxImGW^r^w)] • d fcn 

> / dcoio — hC.c. (52) 

7T . ./n w-S2„fc 
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Note that the two contributions from the Rontgen in- with 
teraction have been collected in a single term as given 
on the second line of the above equation by making use 
of the symmetry (|13[) of the Green tensor. In addition, 
the (vanishing) contributions from the free-space Green 

tensor have been discarded. . . fio r°° 

Next, let us separate the forces F„ into their position- n{ A) — ^ 2-~i J Q 

and velocity- dependent parts. The shifted and broad- k 
ened atomic transition frequencies fl nk are velocity- x ^A^mk • IniG^^r^, 1*^4,0;) ■ d kn ^ ^ 

dependent, so that the first term in Eq. (fS"!?|) also con- lo — Vt nk 

tributes to the velocity-dependent part of the force. 
Again retaining only terms up to linear order in the ve- 
locity, we find 

F n (r A ) + F n (i\4, va), (53) and 



F n (r A ,v A ) = ^ I 

T u JO 



duj LO 



2 n nk (v A )Vd nk ■ ImG (1) (rA,rA,w) • d k 



■ ImG (1) (r A ,r A ,w) • d krt 



l/i£ 
7T 



~ (V'-V)d nfe -v A x [V xlmG^(r Al r Al io)] ■ d kn 

dcoLo hC.c. (55) 



LO — Q 



nk 
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where Vl nk = O nfe (r A ) and n nk (v A ) = n nk (r A , v A ). The f°° ^ u 2 Im (r, r', u) 

velocity-independent force (|5^|) is iust the well-known CP J Q (uj — fl nk ) 2 

force on an atom at rest. We will in the following restrict ^.00 £ 2 (fi 2 — £ 2 )G^(r r' if) 

our attention to the velocity-dependent force f5j} , which = d£ " fc — - — 2 - — - — - 

consists of three terms: The first, generalised Doppler ^° ^ ^«fc) 

term is due to the velocity-dependence of the atomic tran- + ir[oj 2 G^(r,r',u))\ _„ <d(u) nk ). (57) 

sition frequencies; the second, delay term is associated 
with the time interval between emission and reabsorp- 
tion of the electromagnetic field; and the third, Rontgen 
term is due to the coupling of the current density associ- 
ated with the moving atom to the magnetic field. 

In close analogy to the case of an atom at rest, the Substituting these results into Eq. ([55]) . one finds 
force can be separated into its resonant and nonreso- 
nant parts using contour-integral techniques. Writing 
ImG = (G - G*)/(2i), using the property Jl2} of the 
Green tensor, and employing Cauchy's theorem to trans- 
form integrals along the real axis to integrals along the 

positive imaginary axis plus contributions from the poles, Fn( r Aj va) = ^"(r^, v^) + F^(r^, va) (58) 

one can show that 

ujlm G« (r, r', w) f°° , t £ 2 G^ (r, r', it) 



dui = / d£ ■ „ 2 

1 1 LU — \l nk J £ + \l nk 

7rft„ fc G«(ry,ttnfc)e(w„ fc ), (56) with 



~ f V tr{ [ an (v A , iO + c n (vA, -m ■ G (1) ('A, v A , %)} 

- ^ r^t 2 V(v A • V')TT{[ a ' n (it) + a' n (-it)]-G^(r A ,rA,it)} 
Z7r Jo 



and 



F^(r A ,V4) = ii ^Q{io nk )Q. nk {-v A )[^ 2 Vd nfe • G (1) (r A , r A , u) ■ d fe „]^ =n ^ 
fc 

+ i^ ^Q{Cj nk )[uj' 2 V{yA ■ V')d„ fc • G {1) (r A ,r A ,uj) ■ d kn ]l =n ^ 

k 

+ i/x ^e(w„fe)fl„ fc (V'- V)d„ fe -v j4 x [V x G {1 \r A ,r A ,n nk )] ■ d kn + C.c. (60) 



Here. 



a»(w) 



k 



- n: 



dnfcd^ n 



is the polarisability for an atom at rest and 



a n (v A) w) = - 



(61) 



(W + £l„fc) 2 



(62) 



r 



is the correction to this polarisability for a moving atom 
within linear order of the atomic velocity. 

It is instructive to consider the perturbative limit 
fl nk —> uj nk (i.e., 5uj n / k , T n/k — > 0). The resonant force 
can be represented by its zero-order approximation in 
Suj n / k and T n i k which reads 
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F T n (T A ,v A ) = 2no^2@(w nk )[Suj n (vA) - Su; k (v A )][cj 2 Vd nk ■ RcG (1) (r A ,r A ,uj) ■ &kn\' u=IJJnh 
k 

- /io^eK fe )[r„(v A ) + T k {v A )][i J j 2 Vd nk •ImG (1) (r J 4,rA,w) • d fcn ]^ =( ^ fc 

k 

- 2^ E Q{uj nk ) [uj 2 V(v A ■ V')d„ fe • Im G (1) (t a ,t a , uj) ■ d kn ] 

k 

- 2/i ^] 6(w nfc )w rf (V' - V)d nii • va x [V x ImG (1) (r A ,r A , w nfe )] • d fcn (63) 
fe 

The nonresonant velocity-dependent force vanishes to zeroth order in the frequency shifts and decay rates, in contrast 
to the force observed for an atom at rest. The leading nonvanishing contribution is linear in these quantities and it 
reads 



F nr, „ \ 2/i ^ f°° ^ [fa n (VA) ~ 6uj k (v A )](uj 2 n ~ g) ^ A n iX), r . 

h n (r A ,v A ) = > / 7 2 T g2 ^2 vd nfc • G y >(r A ,r A ,i£) ■ d kn 

n k Jo K^kn + S ) 

, 2/Xo ^ Z" 00 2 ^kni^n + ^k){^> k n ~ 3f 2 ) „,> , -,m / . ^ , 

+ y. d ^ 1 2 i ^3 V ( VA ■ V ) d " fe ■ G \ V A,*aA) ■ dfcn 

^ k JO K w kn + ^ ) 

-^E /°°d^ 2 7" 2 (r " + ^ fc) (V- - V)d nk ■ v A x [V x GW(r A ,r A> ig)] ■ d kn . (64) 



For an isotropic atom, these results reduce to 



F T n (r A , v A ) = ^ E 6(^nfc)|d„ fc | 2 [(5cj„(v y i) - Sui k (v A )] [uj 2 V TrRe G (1) (r A , r A , ^)]' u=UJnk 



k 



^E e ( w » fc )l d « & | 2 [ r »( v ^) +rfc ( v ^)][ w2vTrImG(1) (rA,r A , W )]^ =tjiifc 



.3 



k 



^ E 6K fc )|d Ilfc | 2 [c 2 V(v A • V) Tr Im G« (r A , r A , u)] ' u 
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k 



2/io 



E@Kfc)w nfc |d nfc | 2 (V- V)Tr{v A x [VxImG (1) (r A 

)]} (65) 
d fe 



and 



3tt 



F«(rx, v A ) = £ |d„ fc | 2 f d^ e '^^y-^ V Tr G« (r., r., if) 

■ E I**!' d ^ 2 " fc " (r " V fc i ( ^3" 3g2) V(v A • V) Tr G« (r Al r A ,i£) 

^El d ^| 2 /°°d^ 2 7" 2 (F :^ 2 fc) (V- - V)Tr{v A x [V x G W(r A ,r A ,ig)]}. (66) 



J kn 

I 

C. Motion parallel to a planar interface Weyl expansion of the Green tensor 

/d 2 fcii / 
, ' 2 c' k " ,(p ~ p ' ) G(k||, z, z',uj) (67) 

geometries, in oruer lo gam piiysicai msigiiL, we reauriei 

ourselves to the generic quantum friction scenario of an with r = (p, z) can then be used to calculate explicit 
atom moving parallel [v A — vy = (v x ,v y ,0) T ] to a ho- expressions for the terms that contribute to the veloc- 
mogeneous dielectric or metal of permittivity s(uj) whose ity dependent force. The relevant Weyl components 
plane surface defines the (x, y)-plane (see Fig. [T]). The G(kii,z, z', oj) of the Green tensor for z,z' > are given 
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FIG. 1: Atom moving with a velocity v near a planar surface. 



by [Gi = Gij(k||,z, 



G xx 

G X y 
G X z 

G zz 



2k z 
2k* 



k 2 k 2 k 2 

' s p 'P 

rv || ft || 



nj X rvy fx, z rv x r\.y 



P Pfc 2 



1 

2k : 



. e ikz(z+z') r k z k x 



V k 2 



with 



k^ 



hz 



e{u)k z 



Hz 



h z 



e(u>)k z + k lz 



(68) 
(69) 
(70) 
(71) 

(72) 



being the Fresnel reflection coefficients of the surface for 
s- and p-polarised waves [k 2 — uj 2 /c 2 , k\ — e(uj)id 2 /c 2 , 
^Ct)z = ^fl) — ^y]' ^ n e other components of the Green 
tensor can be obtained by using the reciprocity con- 
dition G(r, r',Lu) — G T (r',r,oj), which translates into 
G(k|| , z, z', ui) — G T (— k|| , z', z, uj), and the replacement 
rules Gyy = G xx {k x ^ ky), Gy Z = G xz {k x ^ ky*). 

For the assumed motion parallel to the surface, the 
velocity-dependent shifts and rates vanish Sui^(ta, ~Va) = 
r^(r^,Vyi) = (cf. the remark at the end of Sec. lIII A"| . 
and so do the generalised Doppler contributions to the 
resonant force (|63|) (first two terms) and the non-resonant 
force (|6"6")l (first term). To calculate the delay and 
Rontgen contributions, we require second derivatives of 
the Green tensor as given above. It is useful to note 
that all those derivatives vanish that do not contain an 
even number for each of the cartesian indices (x,y,z). 
For example, terms such as d x d y G X x or d y d y G xz will not 
contribute whereas terms such as d x d y G xy or d v d v G zz 
will. For simplicity, we restrict our attention to the non- 
retarded or near- field limit, where the dominant contri- 
bution to the Green tensor is due to evanescent waves 



with ki z 



ik\\. With this replacement, we have 



r s = 0, 



e(w) + 1 



(73) 



and Eqs. (JHTJ) ((7TJ) lead to 



V(v r V)TrG^(v A ,v A ^) = 



3c 2 v|| e (w) - 1 



16?rw 2 4 e(oj) + 1 ' 

(74) 

(V - V) Tr{ V || x [V x G« (r A , r A , u)] } = 0. (75) 



In the near- field limit, the Rontgen contribution hence 
also vanishes and quantum friction is entirely due to the 
delay effect. 

Substituting Eqs. (J7BJ) and (J75J into Eqs. ([55]) and 
we find the friction force 



K(*A>V\\ 



'II 



4ire z A 



^ 6(w„fc)|d„ fc 

k 

Ime(u;) 



and 



F" r (rA,v,| 



87T£ Z A ^ 



|d„fc| 2 Wfe„(r n + r fe ) 



(76) 



d£ 



J kn 



3£ 2 e(i^) - 1 



(77) 



If we further assume a single-resonance Drude-Lorentz 
model for the permittivity, 



e(«) = 1 + 



2i7tt> ' 



(78) 



with Plasma frequency wp, transverse resonance fre- 
quency lot and line width 7, we find that for a weakly 
absorbing medium (7 <C lop^t) the resonant and nonres- 
onant forces are given by 



F r „(r A ,V||) 



8ne z A 



5 / , 



Q(uj nk )\d nk 
7^(0 



^lk) 



<4) a 



(79) 



and 



F-(rA,v,|) = - 



i2-KEQZ A 



sign(^fc«)(r„ + Tk)uJ± 

UJ S (\ukn\ + ^s) 3 



(80) 



(ujs = \j 'lu 2 ^ + uj'p/2, surface plasmon frequency). 

Let us discuss our results. We first note that in a 
quantum friction scenario of an atom moving parallel to 
a plane surface, a generalised Doppler effect does not con- 
tribute to the velocity-dependent force; this will be dif- 
ferent for an atom moving perpendicularly towards the 
surface. In the near-field limit, the magnetic Rontgen 
coupling becomes becomes negligible as well; it will be- 
come relevant for larger distances. Near-field quantum 
friction forces are hence dominantly caused by a delay 
effect. 

For a ground-state atom, only a nonresonant force 
component (|80|) is present. With both u>ko and T k be- 
ing positive quantities, FJ^r^, vm) is strictly antiparal- 
lel to the velocity and hence presents a genuine friction 
force. Note that this force is proportional to the rates of 



9 



spontaneous decay Tk, the absorption parameters of the 
atom. In the near-field limit, these decay rates are given 

by us, mi, m 



k 



|d„fc| 2 Ime(oj nk ) 
'Qirfie a z\ \e(iv nk ) + 1| 5 



12^£ 4 {uj 2 nk - u*) 



2\2 



(81) 



Inserting this into Eq. (|80[) yields a friction force that 
is extremely short-ranged and falls off as z A s . This 
is in contrast to previous theories [3l that predict a 
2;^ 5 -scaling, resulting from a disregard of the distance- 
dependence of the spontaneous decay rate. 

For an excited atom, resonant forces arise as a conse- 
quence of possible transitions to lower lying atomic en- 
ergy levels. They will dominate the velocity-dependent 
force, in particular if one of them is near-resonant 
with the surface plasmon frequency wg. Depending on 
whether the respective atomic transition frequency io n k 
is smaller or greater than the surface plasmon frequency, 
the velocity-dependent resonant force will either be a de- 
celerating friction force antiparallel to the velocity, or 
it may be an quantum acceleration force parallel to the 
velocity. This can be qualitatively understood from an 
energy consideration. The energy huj n k emitted during a 
downward transition of the atom is resonantly absorbed 
by the surface and leads to the excitation of a surface 
plasmon with energy hcos- The energy difference between 
these two reservoirs leads to a change in the atom's ki- 
netic energy. If the emitted energy is smaller than the 
absorbed one, the atom has to decelerate; if the emit- 
ted energy is greater than the absorbed one, the atom 
will accelerate. Note also that the off-resonant contribu- 
tion (JSOJ) consists of strictly accelerating downward con- 
tributions as well as strictly decelerating upward contri- 
butions which, as before, have a z^ 8 -scaling and can be 
safely neglected with regard to the z^ 5 -scaling of the res- 
onant forces. It is known from previous studies [23| that 
the rate of spontaneous decay increases for atoms in mo- 
tion. This mechanism leads to a more rapid dissipation of 
the internal energy initially stored in the atom, restrict- 
ing the life time of resonant forces; it needs to be taken 
into account in a more quantitative analysis of energy 
conservation. 



IV. EXAMPLES 

In order to illustrate the effect of velocity-dependent 
forces on atoms, we present a selection of numerical ex- 
amples. We will concentrate on forces that are dominated 
by a single atomic transition between a ground state |0) 
and an excited state |1) with frequency to a and (isotropic) 
dipole matrix element d. In this case, the nonresonant 
ground-state force, Eq. (|5D|) . reduces to the simpler ex- 



pression 



F (rA,V||) = - 



32tte z a uj s (uja + ws) 3 ' 



where the decay rate (l8Tj) now reads 



r = 



12nhe z A (u 2 



A 



(82) 



(83) 



The excited-state force is dominated by the resonant 
force component, F^r^Vy) = FY(y a ,V\\) +F\(r A , V\\) 



Fi(i\4,V||) = 



8lT£ Z A 



(84) 



As a first example, we consider a ground-state 87 Rb 
atom moving parallel to a gold surface. We consider the 
lowest electronic transition D2(5 2 Si/2 — > 5 2 P3/2) with 
wavelength Xa — 780 nm [uja = 2.41 x 10 15 rad s" 1 ) 
24j and dipole moment d = 4.23ea = 3.58 x 10 -29 Cm 
25| | . The permittivity of gold may be characterised by 
a plasma frequency ujp = 1.37 x 10 16 rad s _1 and an ab- 
sorption parameter 7 = 4.12 x 10 13 rad s -1 22]. Note that 
the transverse resonance frequency vanishes for metals, 
cjt = 0, so that the surface plasmon resonance is located 
at ojs = u>p/ \[2. With these parameters, we find a deceler- 
ation of the rubidium atom as (ms7Rb — 1-44 x 10 _25 kg) 



= -V||(9.6f 



1 nm 



ZA 



(85) 



The force is extremely short-ranged, and is negligible for 
any reasonable values of the velocity and atom-surface 
distance. 

In contrast, for an excited rubidium atom with the 
same data as above, the deceleration becomes 



V||(5.0 x 10 4 s" 



1 nm 

ZA 



(86) 



In comparison to the ground-state force, excited-state 
quantum friction is strongly enhanced and has a much 
longer range. For an atomic velocity of i> = 200 ms -1 , the 
deceleration at an atom-surface distance za = 10 nm can 
be as large as s = —100ms -2 . Even at za — 100 nm the 
deceleration is still a=— 10 -3 ms -2 . 

Results for other atoms and metallic surfaces can be 
easily obtained by noting that in most cases, the relevant 
atomic transition frequency is much smaller than the sur- 
face plasmon frequency of the metal, hence to a <C t^s ■ Un- 
der this approximation, the excited-state force (f8~4")l and 
the decay rate (|83|) read 
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F i( r ^ v ll)«-^T-^ 



2irsoz A u- 



(87) 



and 



Wi 7 



3irhsoz A to 
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Typical values for the material parameter LUp/j are tab- 
ulated in Ref. [22]. Note that unless the excitation is 
maintained by continuous repumping, the excited-state 
force only acts during a time interval At T" 1 . The rel- 
ative velocity change during this time is approximately 



Av 

v 



Tmv 



3Ti 



2rriAUJAZ 2 A 



(89) 



continuous repumping this force acts only for a very short 
time, leading to a net relative change in velocity 



Av 



S 2 + 7 2 



(94) 



In this limit (w^Cws), the relative change in velocity is 
independent of the strength of the atomic dipole transi- 
tion and all material parameters. 

Upon inspection of the excited-state force (|84"|) one 
notices that this force can be resonantly enhanced if 
an atomic transition matches the frequency of a sur- 
face plasmon resonance. An example of such a close 
match has been pointed out in Refs. [H, [2?| and in- 
volves a sapphire substrate with its principal surface plas- 
mon at As = 12.21 /jm (ws = 1.54 x 10 14 rad s _1 ) and the 
6-D3/2 - ¥ 7Pi/2 transition in 133 Cs with a wavelength of 
\ A = 12.15 /am (uj a = 1.55 x 10 14 rad s _1 ). Near this 
plasmon resonance, the permittivity of sapphire is well 
approximated by [27J 



2i7a-> 



(90) 



with rj = 2.71, uj p = 0.84w s = 1.29 x 10 14 rad s~\ w T = 
0.70w s = 1-08 x 10 14 rad s~\ and 7 = 7.5 x lQ- 3 o; s = 
1.16 x 10 12 rad s _1 ; note that uj s — \/uj\ + 0^77/(7/ + 1). 
With this model, and introducing the atom-plasmon de- 
tuning 5 = uja — k>Si w e find that the force (l76|) and 
the decay rate ([81]) in the vicinity of the surface-plasmon 
resonance may be given as (\8\, 7 -C ojs) 



Fi(rx,vi|) 



V||d 2 n ujp 75 
4ne z A (77 + l) 2 uis (S 2 + 7 2 ) 



2\2 



and 



l2irhs z A (77 + l) 2 uj s S 2 + 7 2 



(91) 



(92) 



With the dipole moment of the abovementioned tran- 
sition being d = 5.85 x 10~ 29 Cm [28|, one finds (mi33 Cs = 
2.21 x 10~ 25 kg) 



a,, = +v,,(7.1 x lO 1 ^- 1 ) 



1 nm 



1 5 



ZA 



(93) 



V. SUMMARY 

We have shown that atoms or molecules in relative 
motion with respect to a dielectric surface experience 
velocity-dependent CP forces. Solving the coupled atom- 
field dynamics for a slowly moving atom, we have found 
an expression for the linearised velocity-dependent force 
on an atom in an arbitrary incoherent internal quantum 
state moving near an arbitrary arrangement of magne- 
toelectric bodies. In general, three effects contribute to 
the velocity-dependent Casimir-Polder force: a gener- 
alised Doppler effect due to the velocity-dependence of 
the atomic transition frequencies, the delay between the 
emission and reabsorption of photons by the atom and 
the Rontgen interaction, i.e., to the coupling of the cur- 
rent density associated with the atomic motion to the 
magnetic field. 

In order to illustrate the general theory, we have stud- 
ied the near-field force on an atom that moves parallel to 
a planar dielectric or metallic surface. Due to the trans- 
lational invariance of the system, the Doppler term does 
not contribute in this case. Furthermore, the delay term 
dominates over the Rontgen term. For a ground-state 
atom the force is a genuine friction force, i.e. a force an- 
tiparallel to its velocity. It is proportional to the atomic 
linewidth and hence very small. In contrast, excited-state 
atoms can be either decelerated or accelerated depending 
on the relative magnitude of their transition frequency 
with respect to the characteristic frequency of the sub- 
strate material. For metals, the force is always decelerat- 
ing while for dielectric substrates with sufficiently small 
surface plasmon frequency, acceleration of excited-state 
atoms can be realised. 

In addition, the force on such atoms is strongly en- 
hanced when atom and substrate are near-resonant. 
Much stronger enhancement can be achieved when the 
atom moves through resonator structures, in close anal- 
ogy to the stationary case [29| . 



Compared with the result (|8B"|) for the excited-state force 
near a metal, we find a significantly enhanced force. Note 
also that, because uja > (i.e., S > 0), the force is 
accelerating rather than decelerating. As a numerical 
example, for a particle velocity of v = 100 ms -1 and an 
atom-surface distance of z A = 100 nm, one would observe 
an acceleration of a = 7 x 10 4 ms~ 2 . As before, without 
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